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ABSTRACT 

3 

The He (n, p)T reaction can be utilized in a proportional counter 
for neutron- energy determination in the region 100 kev to 1 Mev. This 
has already been demonstrated and the reaction cross section in this 
energy region has been measured with a counter. For practical 
application of the counter as a spectrometer, such as in health physics 
work, it is desirable to increase the counter efficiency. A proportional 
counter utilizing an anticoincidence ring to reduce wall effect has been 
constructed for this purpose. Some preliminary tests of the proposed 
system are described^ 
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CHAPTER I 



INTRODUCTION 



The detection of neutrons and especially the measurement of 
their energy have from the first attempts presented a more difficult 
problem than for the other common nuclear particles. This is be- 
cause the usual methods depended on effects resulting from the charges 
of the particles. The experiments by Dee in which he investigated the 
ionization produced in air in a cloud chamber irradiated by fast 
neutrons^ were reported at the same time as Chadwick's announce- 
ment of the discovery of the neutron„ Dee concluded that, if the 
neutron interacts with atomic electrons at all, the process produces 
not more than one ion pair per three meters of the neutron's path. 

If we accept an arbitrary description of neutron energy, E^, as 
intermediate neutrons: 1 kev to 500 kev, fast neutrons: 0 5 Mev to 
10 Mev, then we may say in general that the predominant reaction of 
intermediate neutrons with nuclei is elastic scattering, and that in the 
fast range many other reactions appear, the most important of which 

. [2l 

is inelastic scattering'- , Hydrogen- or methane-filled ionization 
chambers or proportional counters became an important method of 
detecting fast neutrons because the neutron can impart practically all 
of its energy to a proton in a head-on collision, and the recoil proton 
is then the particle causing the action of the counter. However, a 
determination of the neutron energy spectrum by the recoil method 
requires a thin radiator and a double collimation- -fir st of the neutron 
beam and then of the recoil protons--resulting often in too low an 
efficiency. Theoretically, the second collimation, may not be 
necessary because the neutron spectrum can be obtained from the 
recoil spectrum by differentiation, but this leads to large errors since 
we are effectively taking the differences of numbers of about, the same 
size. Nuclear emulsions and organic scintillation counters have be- 
come very important detectors in neutron spectroscopy, especially for 
neutron energies above 1 Mev 
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For energies below 1 Mev it became apparent that exothermic 
nuclear reactions might be useful for neutron spectroscopy. New 
methods were especially needed in health physics applications, where 
any schemes requiring collimation are automatically ruled out. Fast 
neutrons are of major concern to health physicists because of the 
difficulty of shielding against them, and their large relative biological 

[3] 

effectiveness . In 1946 Feld proposed that a proportional counter 

utilizing the {n, p) reaction might be useful for neutron spectroscopy 

in the 10- to 1000-kev range, with the advantage over recoils of 

[4! 

detecting lower energies and making collimation unnece ssary*- . He 
considered the most promising reactions at that time: 



., 14 . .^14 

N (n, p)C , 

35 35 

Cr^ (n, p)S'^'^. 



However, as he expected, these do not very well fulfill the first of 
the general requirements, which might be listed as follows: 

(1) The reaction cross section should be fairly large and its 
variation with neutron energy smooth and free from resonances in 
the energy region of interest (and accurately known). 

(2) The Q value of the reaction should be small, because (E^ + Q) 
is the total energy measured by a counter, and for accuracy 
cannot be too small a fraction of the total. 

(3) There should be no low-lying excited states of the residual 
nucleus to cause an ambiguity in the energy of the reaction particles. 

We can look for these requirements to be met over appreciable 
ranges only in light nuclei, in which the energy levels are more 
widely spaced, avoiding resonances and excited states. Also, for 
small Q values appreciable penetration of the Coulomb barrier for 
(n, p) or {n, a) charged-particle reactions would be limited to light 
nuclei 

Two {n, a) possibilities are 

Li^ (n, a)T, 

B (n, a)Li^. 
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The first has a large Q (4.78 Mev) and a resonance at = 0.Z7 Mev. 
The second, widely used for thermal neutron detection, has an ex- 
cited state and fairly large Q values (2.79 and 2.31 Mev). 

With supplies of helium- 3 becoming more available and a likely 

possibility being He p)T, helium- 3-filled proportional counters 

came under consideration. In 1950 Coon showed that the cross 

[51 

section in the l-Mev region was about I barn , comparable to the 
scattering cross section of hydrogen. The thermal cross section for 
this reaction is large (5400 ± 300 barns)^^^. These types of reaction 
cross sections obey the l/v law for a wide range, but intermediate 
neutrons may depart from the l/v law. 

Some work on helium- 3-filled proportional counters has been 

[7] 

done by Batchelor and others at Harwell, England*- . They used a 
small amount of helium- 3 admixed with xenon or krypton as the 
counting gas and stopping gas to reduce wall effect. Also, a small 
amount of carbon dioxide was added to stabilize the multiplication 
process. In their final filling the pressure of krypton us^d was l64 
cm of mercury. Since they found that, for E^ - 1 Mev, a pressure 
of 8 atmos of krypton is required to reduce the wall effect to 10%, 
they set up a computer program to correct for wall effect. Their 
work remeasured the neutron cross section in the range 120 kev to 
I Mev more accurately than previously, showing a much flatter 
energy dependence than the l/v law. 

They also pointed out a basic source of ambiguity in energy 
determination other than counter defects. As we show in Chaper VIII, 
the maximum recoil energy of a He nucleus is ^3/4 E^. The 

true energy counts run from Q for E^ ” 0 to E^ + Q„ Thus, for 
3/4 E^ = Q, the recoil spectrum begins to overlap the reaction 
spectrum. On the basis of Batchelor's value of Q as 770 kev, this 
means that ambiguities exist if neutrons of greater energy than 1.03 
Mev are present. The scattering cross section up to about 20 Mev is 
comparable in size to the reaction cross section of interest, but of 
course only a fraction of those neutrons scattered deliver a large 
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portion of their energy to the He recoils. In any case, the presence 

of high-energy neutrons reduces the accuracy. 

Gaseous counters have notoriously low efficiency, and as 

helium- 3 has become more plentiful and less expensive, it seems 

worthwhile to investigate the feasibility of increasing the counting 

rate by increasing the quantity of helium- 3 in the counter. As we go 

to higher pressures, special problems arise, such as the demands 

[ 8 ] 

for gas purity, higher voltages, and better insulation . Also, the 
container must be more massive, and proper consideration should 
be given to the effects of the masses upon the spectral distribution 
of the neutrons, primarily caused by inelastic collisions (trans- 
parency correction). For these reasons it would be desirable to 
use helium- 3 alone as the primary filling, utilizing it for the count- 
ing and stopping action. However, we soon find it very difficult to 
decrease wall effect to a reasonable value by helium-3 pressure 
alone because stopping powers decrease with lower atomic numbers. 

For instance, the following ranges of a 5.3-Mev alpha particle at one 

[9l 

atmosphere have recently been measured . 

helium- -Z1 cm, 
krypton--3.0 cm, 
xenon- - Z. Z cm. 

For this reason it is desirable at the same time to determine the 
practicality of reducing the wall effect by surrounding the sensitive 
volume of the counter with an anticoincidence ring. 

Since helium- 3 is quite valuable, a fairly elaborate filling 
system must be devised to insure no loss of gas in filling and sub- 
sequent recovery. The purity requirements are not so great as in 
noble gas scintillation counting, in which a small amount of contam- 
inant such as oxygen or other polyatomic gases may completely 
quench the useful scintillation property, but some purification arrange- 
ments will probably be needed in the filling system to achieve good 
energy resolution and to remove impurities which inevitably get in 
to the gas during filling and recovery. It is important that the 
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loss of helium-3 be kept small during the purification process. Purifi- 
cation procedures for noble gases using metallic calcium have been 
applied successfully for some time in removing all electro.- negative 

[ ] 

impurities, nitrogen, hydrogen, and to a lesser degree carbon dioxide 
Also, the method of physical adsorption of impurities in an activated 
charcoal trap cooled to liquid nitrogen temperatures has been used 
successfully in a helium-3 system, and has the advantage of greater 
ease and convenience^ ^ ^ ^ . However, it appears there may be a special 
problem due to the radioactive background from tritium contamination 
of helium-3, and this problem is increased by using higher helium-3 
pres sure 

It was felt that useful information about these two techniques, (a) 
using helium-3 as the detecting medium and primary counting gas and 
(b) reducing wall effect with an anticoincidence ring, could be obtained 
by testing a counter with a filling of ordinary helium. The natural 
abundance of helium-3 is about 0.00013, so that it is not ordinarily 
noticeable. Owing to the large cross section at thermal energies, it 
can be detected in tank helium by use of a nuclear reactor as a source 
of neutrons, and this has been suggested as a way of surveying world 

r 1 z] 

helium supplies for heiiim-3 content . Because it is a double- 

4 

magic-number nucleus. He is very stable, and the detection of 
neutrons in a helium-filled counter is due to the recoil alpha particles. 
Helium-3 and helium-4 have the same electronic configuration, and 
we can obtain information about the counting characteristics of one 
by experimenting with the other. 
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CHAPTER II 



counter s y s tem 

A method of reducing wall effect in proportional counters is to 
detect the particles that escape from the sensitive volume and exclude 
them from the measurements. To accomplish this it is desirable that 
the walls of the main counter be transparent and completely surrounded 
by transparent- walled counters in anticoincidence. The feasibility of 
Geiger- Mueller counters with transparent walls --so that wall thick- 
ness would not limit the radiation to that above a minimum-- was 
investigated in 1944, and it was found that symmetrically placed wires 

[ 13 ] 

could approximate a cylinder^ \ Single counters of triangular type 
were constructed with good results. However, when two such counters 
were placed close together, there was a deleterious effect due to 
distortion of the field at the central wire. Better results were obtained 
with a regular hexagonal shape. 

COUNTER CONSTRUCTION 

A counter has been constructed in which a layer of gas surround- 
ing the main counter is in effect turned into a separate detecting device 
by a circular array of wires. The geometry is shown in Eigs. 1 and 2. 

The inner circle of wires, grounded to the cylinder, serves as the cathode 
for the main counter and as the inner cathode wall for the ring counter. 

In the outer circle of wires alternate wires serve as cathode and 
counting wires for the ring counter. Since there are 36 wires in each 
circle, we effectively have the main counter surrounded by 18 separate 
counter s. 

An arrangement similar to this, referred to as a "wall-less 

counter", has been applied to low-energy beta spectroscopy and proved 

f 14] 

to be a practical instrument*- ** . As pointed out, this prevents 
escaping high-energy particles from having an effect on the low-energy 
end of the spectrum, but the effective sensitive volume now depends 
on energy. However, inaccuracies because of this circumstance will 
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Fig. 1. Arraagement of counter. 
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be much less than from uncorrected wall effects Another application 
of this type of counter geometry is the so-called '’Giles telescope" 

f 15] 

neutron spectrometer \ In this arrangement the counter is long 
and narrow, and neutrons are collimated along the long axis. The 
requirement of a thin radiator is obviated by using gas recoils and the 
anticoincidence method to prevent registration of pulses that recoil 
oblique to the chamber axis. Monoener getic neutrons then have a 
pulse-size distribution with a well-defined peak. A recent applica- 
tion of this scheme, utilizing hydrogen and propane fillings at 
various pressures, has resulted in energy spreads of 10% for neutron 
energies from a few hundred kev to 10 Mev^^^^ . 

Some of the constructional details of our counter follow. The 
counter cylinder has l/l6-inch-thick brass walls and is 4 inches in 
outside diameter. The total gas volume less the inlet piping is 1.93 
liters. The vessel has been pressure- and leak-tested with helium at 
20 atmospheres, making operating pressures up to 10 atmospheres 
possible within safety requirements. All wires are stainless steel of 
3-mil diameter. The central wire is supported by kovar insulators 
and is electrically available only at the base. Its length between 
thickened supports is 8,5 inches. If we assume that, owing to a fall- 
off in field strength near the ends, only 90% of this length defines the 
sensitive volume, the sensitive volume of the main counter is 888 cc. 
The ring-counter anode wires and all cathode wires are each strung 
from a continuous piece of stainless steel wire and threaded through 
holes in ^-inch-thick lava insulators at top and bottom, the anode 
wire system being electrically available through another kovar at the 
base, A major drawback to this method of stringing is the difficulty 
of repair in case of wire breakage. The volume defined by the ring 
counter is 600 cc. 

The anticoincidence ring also cancels out particles that enter 
the counter from the outside or are emitted owing to contamination of 

the brass walls, Korff says a well-cleaned surface should emit 0,36 

2 . f 17] 

alpha particle per cm per hour within a factor of ten , and 
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Sharpe gives 0.09 alpha particle per cm per hour for commercial 

f 1 8 1 

brass^ . In normal operation the counter will probably not be 
sensitive to single secondary electrons knocked from the walls by 
gamma rays, but unwanted counts due to pile-up will be more likely 
in the ring counter, where the density of secondary electrons will be 
greater. A proportional counter with walls of copper tubing l/3Z 
inch thick has been used for counting gamma rays of energy up to 1.6 



Mev 



[19] 



The gamma-ray attenuation ‘'tenfolding length" is 
2 



approximately 60 g/cm for all materials in the energy region where 
Compton scattering is the most important process This extends 

from less than one to several Mev for all elements. Multiplying by 
density, we find our brass wall is 1.34 g/cm^ thick. The tenfolding 
length converted to a l/e length is Z6.1 g/cm . We have 
exp (- 1 . 34 / 26 . 1 ) = 95 %, therefore about 5% of all gamma rays imping- 
ing on the counter wall will be converted to secondary electrons. Of 
course, many of these electrons will be absorbed in the wall before 
reaching the counter gas, but gamma- ray counting efficiency would 
be reasonably high if the counter were operated at high gas 
multiplication. 

The voltage at which a counter under given conditions is operated 



determines the gas multiplication. For a given counter Rossi and 



Staub have shown the functional relationship of multiplication 



[21] 



to be 



M = M 




Pa 



where 

a = wire radius, 
b = cathode radius, 

P =gas pressure, 

V =»voltage across counter. 

Since multiplication usually takes place in a cylindrical region of 
small diameter near the wire, we might assume our two counter 
systems identical for purposes of figuring multiplication, except for 
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the difference in In b/a. The multiwire potential problem would be 
very difficult to compute, but if we take an average b value for one 
of the ring-counter wires (averaged over measurements to the inner 
circle of wires, to an adjacent wire in same ring, or to the wall), 
for the ring counter we obtain b/a “ 225 and for main counter b/a “ 1000. 
The ratio of the logarithms is 1.27, therefore for the same value 
of gas multiplication, this indicates that the main counter should have 
a higher voltage by a factor of lo27. Of course, the determining 
factor on multiplication of the ring counter will be to make it as sensi- 
tive as possible in order to detect all particles that leave the main 
counter, but not so sensitive that a large number of background counts 
keep the desired counts from the main counter gated out. 

ELECTRONIC CIRCUITRY 

A block diagram of the electronic circuitry used with the 
counter is given in Fig. 3. The blocks below the dotted line complete 
the neutron- spectrometer system, but investigation of the operation 
of the system was made without this by scaling both coincidences and 
anticoincidences, A list of units actually used is given in Table I. 

One channel of the dual- channel variable gate and delay unit can be 
used. for the variable-gate block, and the use of this unit for the 
pulse - shaping function could be replaced by simpler electronic cir- 
cuitry. This unit is very useful, however, for experimenting with 
different values of discriminator setting and of delay and length of 
output pulses. If two pulses arrive at the pulse- shaping block at 
about the same time, we want to insure that there is no anticoincidence 
pulse to gate the pulse-height analyzer on. To do this we delay 
channel No. 1 slightly so that the channel No. 2 output pulse will 
definitely be at the coincidence unit before the channel No. 1 pulse 
arrives.- Also, channel No. 2 pulse is made longer so that the 
coincidence unit will see it during the entire time that it sees the 
channel No. 1 pulse. The count on Scaler 1+2 indicates the number 



- 11 - 



I ^ «• 

m — W9 (P " *^ M# iw# 1 1 » )l «»ii<‘i> !••* 

«Am #*o ••• -> » * «iik' -«*- >«k««ii» ffom^^twam 

«*eifc_i«fc l’ ^ - - ' «»“* ti 

41 "'*n-* *• «JUOr^ 

II » iVii • « |i>^ «XP*» 4 ** •" <)(• 

-«*■»■ iS’- >»| 4»| |4>iP*'4 ^Hw.. <o/ ' 

IK* iV- MBifl lifTt»« tlMW- •« «■ iMrifll 

|iil**-V^ 1* •** • 

■ » 1^1 1— II •■» ,r ' I •* • I » I M '««• <■» I ■»««• Ok «Bu> i«i (mr^ 

(^koriiv-oK mFs^ •t — .. .«*^ ««••-» * T _!•(■•»• 

- o4ii-iiiu»- * - 4M 

MC^ *1 Bh 4W firmiNl *f>^ i ^ii # mm'V >• ' «4 i» 

I itiMf 4^ <«M|ii 

*k . I mt. I ._ -ibU-Bfc • I ut ^ I 

» « W J J!f^m I i^%if Jf^ 

• immt - -lISTr !>• «^Mti I ^I'l 

«• ««nr«i l« « f I f i«i| « <|m(H 

^ ^ oJuUlli^l^ AWtMJk. y- - ^ «M»| ai, yn» 

'Oii «w*iii«iii •*%•»• I •> H^|H4pNpw"i 
^o«ifiiO» o»i* f >r«n |«iwia» -i»| >i«l • 4«l|iijn4]r <flpfi»i«aaM«i 

P" • — *«r^ . 

«|T«4 l > v * > * ’>' ; • • in44 

!»’• ) li' • »aO # • II '"'<1 iTlHBSfc <»- o «n -^juw» 

•*. f I fJLi^v •;! 4a»r I 

I 



I 



I 





MU-15.210 



Fig. 3. Block diagram of counter electronic circuitry. 
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Table I. Units used in counter electronic circuitry. 



F unction 

Regulated High Voltage Supply 



Meter and Divider 


Panel 


Linear Amplifier, 


ch. 


# 1 


Preamplifier, ch. 


# 1 




Linear Amplifier, 


ch. 


# 2 


Preamplifier, ch. 


# 2 





Pulse Shaping 
Coincidence Unit 
Scaler 1 + Z 
Scaler 1, Z 



Unit 

Northern Scientific Co. , ser. # 1 

UCRJL 3 kv Divider Panel 

UCRJL Mod 5, Dwg. 1X5334 

UCRL Dwg. IX4353D 

UCRL Mod. 3, Dwg. ZT4404 

UCRL Dwg. 3TZ393 

UCRL Variable Gate and Delay Unit 
Dwg. ZT8084F and ZT8154 

UCRL 8 Channel Quadruple Mixer 
Dwg„ 3T5894-1 

UCRL 10Z4 Scaler, Mod. Z 
Dwg. 3T8934 

UCRL Z Channel 10Z4 Scaler, Mod. 1 
Dwg. 4T6895 
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of pulses from channel No. 1 not registered by the pulse-height 
analyzer. The versatile variable gate and delay unit, used with one 
input pulse and in conjunction with the coincidence unit, can also be 
used as a single-channel pulse-height analyzer if desired. 

IONIC MOTION 

With this size of counter one may have to worry about the time 
lag between formation of the pulse by the main counter and by the ring 
counter when ionization occurs close to the edge of the main counter. 
We can compute the travel time for ions if we know the mobility and 
how drift velocity depends on mobility and the field. Sharpe has 
developed convenient equations for computing ionic motion in 

r 1 8] 

chambers^ ^ . In general, in cylindrical or spherical fields, the field 
is 

Vr'^ 

€ ^ .^P — — , 

In y 

where 

y = b/a 
r - radius, 

V “ voltage across counter, 
and for a cylinder n 1. 

An ion crosses the system from one electrode to the other in a time 
given by 

T = / [(np+U K (V/ln y)P] , 

where p - ^ for electrons and 1 for ions in the equation relating drift 
velocity w to mobility K, 

w = KfiP 

We are concerned with the time for an electron to travel across the 
counter because most of the pulse is formed near the wire owing to the 
concentration of the field at the wire whether there is multiplication or 
operation as an ion chamber. If we take a case of V - 1000 volts, 

P = 2.5 atmospheres of helium, and use Sharpe*s value of 
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K = 3.5 • 10^ cm^' \l volt - sec.-^Jatmos , we find for electrons 
T =7.55 |JLsec. However, we expect that the mobility of a helium- 
carbon dioxide mixture may easily be greater by a factor of ten, 
based on the way CO^ increases the drift velocity of electrons in other 
noble gases. This is explained for argon on the basis of the Ramsauer 
effect and the high first excitation potential for argon, 11.5 volts, in 
comparison with the low excitation levels in a molecule ^ . The 

first excitation potential of helium is even higher, 19.77 volts, 

Extensive work has been done on collection times in argon, krypton, 
and xenon, and it is reported that once a fraction of one percent of 
CO^ is added, the effect of original polyatomic impurities in a noble 

f zzl 

gas may be assumed insignificant^ ^ . The transit time computed 
above is reduced by the same factor as that by which the mobility is 
increased, and should not be too long for the coincidence requirements. 

The transit time for a He^ ion in pure helium to cross the system 
from wire to wall, computed for the same conditions as for the electron 
and on the basis of a mobility of 10.8 cm. ^/volt- sec-atm. is 1850 psec. 
This time is not too important, since the pulse can be shortened by 
shortening the amplifier time constant. 
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CHAPTER III 



ENERGY RESOLUTION 

A discussion of energy resolution depends on the complete theory 
of proportional- counter action. This has been well covered in several 
books and articles^ ^ . The attempt here is 

to point out important factors concerning the energy resolution of this 
particular counter, using developments from any of the basic sources 
but modifying the terminology as necessary in order to be consistent. 

INHERENT LIMITATIONS 

Inherent fluctuations in a proportional counter are from (a) the 
number of ion pairs released by monokinetic radiation, and (b) the 
size of the avalanche each ionization electron produces. From 
theoretical work based on the assumption that the first was a Poisson 
distribution (it was found experimentally to be smaller than this) and 
expanded upon experimentally, it was shown that relative variance of 
total size of avalanche is approximately given by 

where 

Vp = variance (mean square deviation) of size of avalanche, 
mQ = mean of number of ions initially liberated, 
m^ = mean of number of ions produced in avalanche by one 
electron. 

We note that this inherent spread is no longer a significant 
factor at the energies at which we will be working because in the 1-Mev 
energy range and at an assumed mean energy to produce an ion pair 

4 

of 30 ev, mQ is greater than 10 . These theories were useful, and 
experiments were done down in the ke\r region where mQ was a smaller 
number. For our magnitude of ^tIq, the relative variance would be 
around 0.01%. Since experience has shown that the relative spread of 
pulse height (relative standard deviation of pulse height from a mono- 
kinetic radiation) is at least 1% or Z% (owing to mechanical defects of 
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a counter or defects in gas filling) at these energies, the inherent 
spread is not a significant factor. The mean energy per ion pair, 
which determines niQ, is expected to be practically constant over the 
(E^ + Q) range. This may not be true in very pure helium, in which 
the mean energy has been measured as 41.3 ev, in contrast to Z9«7 ev 
when a small fraction of another gas is present. A possible explana- 
tion on the basis of charge exchange and the metastable states of helium 

r 18 1 

IS given by Sharpe with reference to experimental findings^ ^ . 



NEGATIVE-ION FORMATION 



An important gas defect is negative-ion formation due to electron 
attachment to a gas molecule. It does not occur in pure noble gases 
nor in N^, CO, CO^, or CH^. However, it is likely to occur in 

the collisions of electrons with the halogens, O^, or water vapor. In 
the development of Eq, (III-l) the following equation was derived for 
the case where a fraction h® of ionization electrons reaches the wire 
without forming negative ions: 



I V 

, vp , . + J_ 

(h'mAmo)^ ^2 





(III- 2) 



where = variance of number of ions initially liberated. For h* - 1 
(no negative-ion formation), this equation leads to the previous 
approximation that relative variance is about l/mQ when we use an 
experimental result that - m^/s. That is, the relative variance is 
approximately 



1 



m 



0 



j. 

3 





(A bibliographical review on the development of Eqs. (III-l) and (III-2) 
is given in Reference 26). 

Thus, for accuracy in low-energy work it is recommended that 
the gas be purified to remove negative-ion-forming impurities. How- 
ever, since the factor in brackets in the expression for relative 
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variance IS multiplied by the reciprocal of we expect that in our 

case the purification requirements may not be so stringent. This can 

[is] 

be roughly checked by an approximate expression^ ^ for the con- 
centration C of electronegative gas in a mixture in terms of the gas 
pressure P (in atmospheres), and Shubweg 6 (in cm) based on gas- 
discharge theory: 




(III- 3) 



The Shubweg 6 is the product of the mean free life of the electron and 

its drift velocity under the applied field. Values r epre seui tative of 

the case in many gases were assumed for the electron mobility and 

probability of attachment per collision, h, for the gas. 

Suppose we wanted a gas purity that would make h* in our 

expression for relative variance equal 0.99. If an electron started a 

distance r from the central wire of a counter, the probability of its 

reaching the wire without attachment, by the definition of the 

^r/g -r/6 

Shubweg 6, would be e ^ . Therefore we set e ^ ~ 0.99 or 

6 ^ 100 r. For a counter in which 3 cm would be a representative 

value for r and at a pressure of Z.5 atoms we get, from Eq. (Ill- 3), 

C = 0.01%, a fairly stringent requirement. If we try h* = 0,9 and 

make the same assumptions, we get C ^ 0.1%, or 99.9% purity 

required. Now, going back to the approximate expression for relative 

4 

variance and assuming m^ - 10 , we can solve for a value of h‘ that 
would give a relative variance of, say, 0.1%. This turns out to be a 
seemingly very easy criterion, h® 0.179. However, making the 
same assumptions as before, we find C is about 1%. But in the 
equation for C the assumptions for electron mobility and h, which 
exhibits resonances for some gases dependfent on electron energy, 
could both conceivably be off by almost a factor of ten, and considering 
also other approximations, we could easily be back to our 99.99% 
purity requirement. The parameters in the gas-discharge theory are 
not accurately known for all the pure gases and cannot be predicted 
for mixtures, especially if the percentages are not very accurately 
known. 
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RADIOACTIVE BACKGROUND 



The problem is different if there is an appreciable tritium 
impurity in the helium- 3 supply. The loss of energy resolution is then 
due to the pile-up of low-energy background counts due to the radio- 
active decay of the tritium. This decay has an end-point beta energy 
of about 18.9 kev 

To estimate the required purity of helium- 3 with respect to 

[2 1 1 

tritium we must first analyze the effect of pile-up in generaU 
Consider the piling up of square pulses of uniform height W and equal 
width T . A pulse height n W is produced if n pulses occur within a 
time T. If a scaler discriminatpr level were set at n W, an erroneous 
count would then be registered. We can get approximate values for a 
counter by taking r equal to the equipment resolving time. 

If the occurrence rate of events in the counter each giving rise 
to an energy E has a mean rate of the mean number of events 

during a resolving time is n^r and the standard deviation is for 

n^T » 1. Thus, the energy dissipated per resolving time is 
n^r E±N/nQTE . Now suppose the threshold energy of the system is 
such that it takes an energy of nrE to produce a spurious count. Then 
the following expression gives the number of standard deviations by 
which the number of counts during a resolution time must exceed the 
mean number in order for a spurious count to be produced: 

(tiT - ^inQT - [(n - Qq) / n/'V. (III-4) 

From this we can estimate the probability p^ of receiving a spurious 
count during the resolving time for a given and t. The spurious 
counting rate is given by I - p /r. If we estimate t as 0.5 fisec and 

wish to keep I less than 0.1 count per minute, p must be less than 

/ -9 ^ . . 

IsT - 5/6 • 10 . This probability corresponds to a deviation above 

the mean between 5 and 6 standard deviations. Let us require that 

Uq be less than n by 6 standard deviations. The average energy of the 

tritium beta spectrum is about 5.7 kev, therefore let E equal this 
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value. The range of beta particles of this energy is such that most 
of the beta energy will be spent in the counter gas, as shown in 

3 

Chapter VI. Since the Q of the He (n, p) T reaction is 770 kev, a 
value of nrE less than this will not be detected as a spurious count. 

Of course, a value of nrE that is an appreciable fraction of 770 kev 
destroys the accuracy of neutron-energy determination if enough of 
the pile-ups occur in coincidence with valid neutron counts. However, 
the number of accidental coincidences between two counting rates is 
proportional to the product of the rates and the resolving time; in 
this case the true and spurious rates are expected to be small enough 
that the problem of modification of true pulse-height counts will not 
be as serious as the number of spurious counts recorded compared to 



the number of true counts. If we set nrE ~ 770 kev, n is approximately 
3 • 10^. Using Eq. (III-4), we obtain 



6 ^ 



10® - n. 



T7T 



0 



(0.5 • 10 



- 6 , 1/2 



^8 



Solving, we have n^ ^ 1.86 * 10 disintegrations per second. 

If this is the allowable mean disintegration rate, we can solve 
for the allowable number of tritium atoms in the counter from the 
disintegration equation: 



dN/dt = - (0.693/ ^ - Hq , 

where 

N ~ number of tritium atoms 

Yjz ” half life == 1Z.Z6 years. 

17 

This gives N - 1.03 * 10 tritium atoms. Or if each molecule is 

1 6 

diatomic, the allowable number is 5.65 • 10 ' molecules. The 

sensitive volume of the central counter is V ” 888 cc, and at P - Z.5 

^ ZZ 

atmospheres, the total number ofmnolecules is D P V - 5.66 • 10 

o ^ 19 

molecules, where D ^ Loschmidt*s number at 15 C, or Z.55 * 10 
This gives a tritium composition of about 0.0001%. This looks like a 
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very serious purification problem if there is significant tritium contam- 
ination in the helium-3 supply. 

SATURATION EFFECTS 

It is apparent that for high-energy particles and high enough 
multiplication, saturation effects near the wire should eventually 
destroy the proportionality of the multiplication process. Experiment- 
ally, the ratio of pulse sizes has been found accurate to 1.5% for M 

r 28 1 

less than a critical value of M^(E) at which saturation effects begin^ 

For the counter used the experimenters found 

E • M^(E) = constant = 10^ ev, (III-5) 

Using this for our counter application, where the maximum (E^ + Q) is 
1.77 Mev, and allowing a safety margin of a factor of two, we have 

M (E = 1.77 Mev) ^ 28 . 
c 

If we also check expected charge-to-length ratio in avalanche compared 
to charge-to-length ratio on the wire, a ratio of 4% corresponding to 
Eq (III- 5), it appears that adoption of the equation for this counter 

- •'W 

should ensure no significant loss of proportionality due to saturation. 
High values of multiplication are also not desired, owing to increased 
voltage - stability requirements, because the curve of multiplication 
versus counter voltage in most cases starts out rather flat but soon 
rises into a steep ascent. 

POSITION OF IONIZATION 

Because of the long travel time of positive ions to the cathode, 
amplifier time constants must be used that ignore part of their con- 
tribution to the pulse. This makes the size of the pulse dependent on 
the position of ionization. A general rule is to use M > 10 to avoid 
this effect when counter radius is 100 times wire radius. Let us 
investigate the effect for our counter, for y = b/a - 1000, and for no 
multiplication. On the basis of the expressions for fraction of total 
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induced voltage due to electrons and positive ions, Staub shows that 
the fraction of the volume of the chamber in which the ionization can 
originate so that more than half of the pulse is caused by the fast 
electron collection is b/ (b+a), since the contribution to the pulse is 
evenly divided when r == n/ (a • b) is the position of original ioniza- 
tion Or this is y/(y^I)s which equals lOOO/lOOl, which seems to 

indicate that multiplication might not be necessary for large y values. 
This assumed position corresponds to r - 21. 6a for our counter. As 
a further check we assumed the ionization took place at r - 316 a, 
and calculated the fraction of the pulse due to electron motion 
(ionization in about 90% of the volume would give larger pulses), 
using the assumptions about ionic motion used in Chapter II and 
some other equations from Reference 18. This fraction was 0.835. 
Amplifier time constants in reality would be long enough to register 
some of the positive-ion motion, which would improve this fraction. 
However, starting from r - 3l6a under the same conditions and 
assumptions, we calculate the ionic motion to contribute only about 
0.002 to the total pulse in a collection time of 5 psec. Thus we con- 
clude some multiplication is necessary to remove effects due to 
position of ionization. 



END EFFECTS 

With no preventive measures there would be an abnormally 
high electric field at the ends of a counter wire. The loss in 
resolution thereby can be partially erased by having the wire 
shielded or thickened near the ends to give a low field. The variation 
of M along a counter wire has been measured and found to be import- 
ant at 4 distance on the order of a counter radius from the thickened 
[291 

ends . In a high-pressure proportional counter designed to be 
used for neutron spectroscopy it was found that the pulse-height 
distribution varied radically over the extent of the counter in a pre- 
liminary design without field tubes^^^^. The loss of resolution com- 
pared with the resolution at the center, was significant when the 
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investigators naoved a narrow pencil of thermal neutrons with which 
they were probing the counter 0.5 inch away from the center. The 
wire was about 8 inches long. The simplest way of overcoming end 
effects is to use a long counter - -that is, long compared with the 
radius. Of course, this is not always practicable. Our counter, 
with wire length Z1.6 cm and counter diameter 7.62 cm, is not a 
very long counter, and we might expect considerable loss of resolu- 
tion due to end effects. The end effects can be eliminated by mechan- 
ical methods, which add complexity to the counter. A divided counter 
was constructed in which a bead divided the wire into two unequal 
lengths^^^^ . Subtracting one spectrum from the other eliminates 
end effects. Field tubes extending over a guard tube eliminate end 
effects when field-tube potential is adjusted to a proper value and 

field and guard tubes are correctly located. Field tubes were proposed 

[3i] [321 

in 1951 and used the same year*- ^ to investigate the low-energy 
spectrum of tritium down to approximately 200 ev. 

WALL EFFECTS 

Wall effects as outlined for a gas recoil chamber by Rossi and 

[2 ll 

Staub are pertinent for charged-particle reactions, since in both 
cases the loss in resolution results when the particle that causes the 
ionization in the gas does not expend its full range in the sensitive 
volume. Three cases are: (a) The reaction takes place in the 

sensitive volume but one of the charged particles produced crosses 
the end boundary; (b) The reaction takes place in the sensitive 
volume but the charged particle hits the lateral wall; (c) The 
reaction takes place outside the sensitive volume, but one of the 
charged particles enters the sensitive volume. The anticoincidence 
ring does not help in those instances in which the charged particle 
crosses the ends. Thus, it is desirable to have a small ratio of area 
of ends to lateral area, or- -what is the same thing--aTong counter, 
as in reducing end effect. Pulses reduced in size by the wall effect 
contribute a continuum of energies from zero to (E^ + Q) to the 
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pulse-height distribution. Without the use of an anticoincidence ring or 
a magnetic field to curve the paths of the charged particles, the wall 
effect can be reduced only by increased size of the counter and in- 
creased gas pressure. 



OTHER EFFECTS 

Inaccuracies in energy resolution can also come from variation 
in applied voltage, variation in position and diameter of the wire, and 
dust particles on the wire. These are discussed in the general 
references. Any alone could be very detrimental, but it appears 
that all can be avoided with careful techniques. 
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CHAPTER IV 



GAS FILLINGS 

The counter was first tested by filling from a bottle of argon 
plus 4% carbon dioxide. The filling and evacuation system used had 
an oil diffusion pump and liquid nitrogen trap, but no special purifica- 
tion procedure or trap was used in the gas-filling line. Satisfactory 
operation was obtained in counting gamma rays at pressures of one, 
three, and five atmospheres. Since the secondary electrons have a 
continuous distribution of pulse heights, we would not expect a flat 
plateau in the curve of counting rate versus voltage in the proportional 
region, but we do desire a slope that does not place too exacting 
requirements on voltage regulation. At 5 atmos we found about 
increase in counting rate per volt in the 2600-to- 2700- volt range on 
the main counter. This was at a high value of M; the range of pulse 
heights available was several thousand times the minimum detectable 
pulse height. At 3 atmos in the same voltage range, the change in 
counting rate per volt was about 3%. Of course, the M was greater 
owing to the lower pressure. The ring counter at 5 atmos had about 1% 
increase per volt in the 2 100- to- 2200-volt range. Its counting rate 
was higher owing to its proximity to the wall. 

The helium filling was from a bottle of 99.99% pure helium; the 
desired percentage of carbon dioxide was previously introduced through 
a separate filling line. Both filling lines were run through dry ice- 
acetone slush traps when CO^ was used or a liquid nitrogen trap if 
helium alone was used. {CO^ would be frozen out by a liquid nitrogen 
trap. ) The use of the traps did not seem to have a large effect. The 
CO^ was found very necessary to avoid spurious discharges in the 
counter . 

The stabilizing effect of polyatomic gases mixed with noble gases 
seems to depend on their ability to absorb photons, in contrast to the 
small-photon capture cross section of the noble gases (many of these 
photons are produced by metastable atoms), and also possibly on their 
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ability to receive energy from metastable states that have high 
excitation energies. Colli has investigated ultraviolet photons in the 
decay of metastable argon atoms and found them to be responsible, 
through a photocathodic process, for the starting and maintenance of 

[331 

corona current in cylindrical argon counters*- . The photons were 

produced in the Townsend avalanche on the wire. Another effect of 

metastable atoms is the release of secondary electrons from metal 

[ 34 ] 

surfaces , The positive ions can also eject secondary electrons 

from metals, the requirement being that the energy of ionization be 

greater than twice the work function of the metal, which is about 5 ev 

for common cathode walls. Helium is expected to be more troublesome 

than other gases because of its high ionization potential, about 24.5 ev, 

and the high energies of its metastable atoms, about 20 ev. Evidence 

for the difficulty of interpreting experiments with helium has been 
[35] 

presented^ on the basis of the existence of the helium molecule-ion 
and metastable molecule. It had been shown the He^ ion is present 
at very low pressures, and the concentration increases with pressure 
relative to the concentration of He^. 

The stabilizing effect of was demonstrated in the counter by 

determining the approximate maximum voltage on the main counter 
that could be applied without causing a rapid breakdown when a radio- 
active source was placed a given distance from the counter. Total 
pressure was 5 atmos. The results are listed in Table II. The break- 
down caused by background radiation alone occurs with pure helium 
at 800 volts. 

However, adding more CO^ was not a solution because the 
pulse heights were thereby drastically reduced. It had been found 
previously that in mixture of argon plus 10% CO^ there was a loss of 
pulse height above pressures of 3 atmos, probably due to columnar 
recombination, which is high in CO^. This refers to recombination 
along the track of a densely ionizing particle. It was not possible to 
attain a large enough range of pulse heights above a minimum dis- 
criminator setting for counting recoils from a Pu-Be (a, n) neutron 
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Table II Percentage of in helium versus voltage 

for breakdown on main counter. 



Percentage of 


Critical Voltage 


4 


1500 


6 


2000 


12 


2500 


20 


2800 


25 


3000 


33- l/3 


— 



source without going to a counter voltage at which spurious discharges 
would occur. The situation improved with decreasing pressure. 

A final filling at Z.5 atmos with 2% CO^ was chosen for further 
exper imentation„ The curve of counting rate versus voltage for count- 
ing neutron recoils from the Pu-Be source showed a slope of 0.84% 
change in counting rate per volt in the range of 1100 to 1200 volts on the 
main counter wire. For the outer wire there was a 4.3% change in 
counting rate per volt in the range from 800 to 1000 volts. The 
voltage - stability requirements at 5 atmos at a voltage at which the 
spurious discharges did not occur were not significantly different. 

One advantage of helium as the counting gas is readily apparent-- 
its reduced sensitivity to background gamma radiation. The range of 
the secondary electrons in g/cm^ is about the same, but helium has 
about l/ 10 the density of argon. The low atomic number of helium is 
also an advantage in that the interaction of gamma rays directly with 
the gas is small. 
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CHAPTER V 



COUNTING RATE AND EFFICIENCY 

For an isotropic neutron flux in a proportional counter the count- 
ing rate naay be expressed^^^ as 

I = NV^ f <,(£__) 

where 

N = density of reacting nuclei (nuclei /cm ), 

3 

volume of absorbing material (cm ), 
a (E^) - neutron cross section of absorbing nuclei 

(cm^)j 

(j>^ - neutron flux per energy interval (neutrons/cm - 

sec-erg), 

F(E^| - detection efficiency defined as ratio of count- 
ing rate to rate of reacting events. 

Ideally, for an exoergic reaction and integral detection, F(E^) is a 
step function: zero for < B and unity for E^ > B, where B is 
threshold bias energy for equipment. This step function is modified 
by a geometric factor, G(E^), for wall effect as shown below. For 
differential detection, F(E^) is ideally a difference between two step 
functions with thresholds Bj and B^, For a recoil- detection process 
in which, with no collimation, a neutron of any energy may give 
pulses from zero to a maximum, we must evaluate F(B/E^), which 
gives the fractional number of recoils giving pulses larger than B. 
Thus, for a recoil detector, one has 

F(E^) = G(E^) F(B/E^). 

The term above is a correction factor for perturbation of 

flux caused by insertion of the absorber. This factor is neglected in 
the rest of the discussion. 

A solution of the equation for counting rate would in most cases 

be extremely complicated. Some insight can be gained by considering 

the simpler case of a monoergic neutron flux. Then we can eliminate 

the differential and replace <j)^ with (E^) in neutrons/cm -sec. 

n 
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For neutrons of a given energy this leads to 

I rr N V a(E ) <b(E ) F (E ) . 
c n ^ n ' n' 

The efficiency of a counter is defined as l/q, where q is the 
number of neutrons per second traversing the counter volume; 
q = cj){E^)A, where A is the effective exposed area of the counter. Then 
counter efficiency is 

NVca(Et^) <j)(E^)F(E^) N a (E^) <j>(E^) 

A A • 

If we call V /a = d, the average distance traveled by a neutron in 

^ [ 24l 

crossing the counter*- , counter efficiency is 

e^ = N a (E ) d F(E ] = e^ F(E ) , 

C ' n' ' n^ R ' n' ' 

where e^^ - N O’ (E^) d is the efficiency of radiation defined as the 
ratio of the rate of reacting events to rate of neutrons traversing the 
counter volumeo 

The quantity N is readily calculated by expressing it as D P, 
where D is Loschmidt’s Number and P is the pressure in atmos. 

For the counter, we have d = 1.59 cm if we take the total sur- 
face area of the main counter as the effective exposed area for an 
isotropic flux. K we assume 2.5 atmos of helium-3 and a neutron 
energy for which the (n, p) cross section is one barn, we have 
eR ^ 0.01%. For the counter filled with ordinary helium at the same 
pressure and for a neutron energy for which the scattering cross ^ 
section is 2 barns, we have e^^ ^ 0.02% for recoils. 

The opposite extreme to an isotropic flux is a collimated neutron 
beam. Here the general equation for counting rate is 

I = A <t>E (1 - e'^ '^) F(E^) dE^ , 

n 

where 

A = area perpendicular to beam, 

T = thickness in cm of material parallel to beam. 
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Assuming a monoargic neutron beam as before, we can write 

I = A4.(E__) |1 - e'” '' 1 F (E^) , 

e^=[l-e-N<'(E„)T] F ,E^) , F(E„) . 

and for N a (E ) T « 1, e^ ^ N a (E ) T. 

n tv n 

In this case we see that the radiation efficiency of a cylindrical counter 
in a neutron beam is dependent on its orientation with respect to the 
beam. If the beam direction is perpendicular to the long axis of a 
counter of radius b, solving for the averp,ge thickness of path gives 
( tt/2) b or 5.98 cm for the main counter. Thus, the counter has a 
higher radiation efficiency under these conditions by the factor 
T/d = 3.76. Of course, e^^ would be even higher for a beam down 
the long axis of a cylindrical counter. 

EXPERIMENTAL SETUP 

However, in experimental setups the counter usually is neither 
in an isotropic flux nor exactly in a beam. Consider placing a neutron 
source at point P, a distance c from the counter axis and opposite the 
center of its length, as shown in Fig. 4. If c is not too small and the 
counter dimensions not too large, it seems we could approximate the 
expected counting rate by assigning an average neutron flux to the 
counting volume: 

where = rate of emission of neutrons by the source and 
r^^^ = (c d) , where d = Vb^ + f^ + c^ . Then, 

I = N V (7 Q/ 4ttc d 

(assuming F(E^) = 1 for the time being). The rate for neutrons 
traversing the counter in neutrons per sec^Dnd can then be expressed 
as this average flux times the area of the central rectangle, or 
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Fig. 4. Geometrical arrangement of neutron source and 
counter . 
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This expression was set down to show that this method is equivalent 
to a common method of computing the solid angle subtended by a finite 
rectangular counter^ . For one of the small rectangles of area b f 
in the geometry shown the equation for solid angle given is 

tan = b f/c d 

For a small angle we have tan Q and for the entire rectangle the solid 
angle would be 

4 b f/c d. 

Thus, we have 



q -^4^ / 4iTj- Qg~(Qg/u)- (bf/cd) , 

as above. Also one should note that we have assigned the same 
radiation efficiency in this geometry as for a beam, as we can show 
= l/q = N cr V^/4 bf (still neglecting F(E^). 

But we have V^/4 b f - irb^ • 2 f/4b f = (Tr/2)b, the average thickness 
found for a beam through the counter. 

An arrangement like this was used with a Pu-Be neutron source 
of Qg = 1.54 * 10^ neutrons/sec and d = 21.2 cm. Then we compute 
I/F(E^) = 1662 cpm for a pressure of 2.5 atmos of helium, and using 
an average scattering cross section of 2 barns corresponding to an 
average neutron energy of about 4 Mev. (Actually the Pu-Be neutron 
spectrum is complex from zero to about 11 Mev, and no attempt has 
been made to compute an accurate average value. ) 

If we obtain curves at different values of high voltage, as shown 
in Fig. 5, we can estimate the point where they tend to zero pulse 
height. At this point F(E^) is unity and the experimental value of 
counting rate should correspond to 1662 cpm. It appears to be only 
about 400 cpm, not a close correspondence, but there were so many 
inexactly known quantities and dubious averaging processes that it is 



- 32 - 



-Ill f — *.!«%. mwAiM* t H 

|| «Mm ki>JM»» tll>» liJi 

i ^ t ^ It *Y«^ 

mrnmm- m %tmi* ** ««*• M«m|«ll JU» MSif Ol'X* I **>^ 



|M'#M 





««* f 



ftm « e4 m <n p ii;? 

, li |«i9«rt'ae 



<«■»•«•• ^ dl '^«I ^ 9 ^^ 




ail iai«» <^90i ♦ Ml 
»l| ^ 

<«■;»« itMii ••»• ««•••■» 4 • 



• - lAtia 

y> t — .immu»ui i 

ji*Dujr ‘**1 9i*rt4 «w«i «>tar» f» 

^1^1 'aiifi in m «Afiw !»• »► • ^ 41^# 

«>iaw «|^l J 

4»a ■ M» waif a» 4Mir «4P^Wia> *t» • 

« 9# ^ #;i ^ f • ' # ai 

[^ •■Mil linn ^ «» IK » niM |l•'••f P• irKK»^«Vlll 




COUNTING RATE (cpm) 



1000 - 
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Fig. 5. Integral bias curves for main counter at various 
counter voltages, 

A. 1 200 volts, 

B. 1 16 5 volts, 

C . 1135 volts, 

and using neutron source in arrangement of Fig. 4. 
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hard to draw definite conclusions. The scattering cross section is 
rapidly varying in this energy range and the manner in which the inte- 
gral bias curve approaches zero is not known. Also consider the count- 
ing rate of 175 cpm for the highest counter voltage used at ^ pulse 
height of 10 volts, the lowest discriminator setting for which counts 
were recorded, K we let G(E^) = 0.83--a value computed for an 
average alpha recoil in Chapter VI then, for a 10-volt bias, we have 
F(B/E^) = 0.127, which would indicate a detection of a small fraction 
of the recoil alpha particles. The curves for the ring counter. Fig. 6, 
seem to tend to a much higher zero-point counting rate, but it is 
believed that other events than the recoils are being counted, since the 
volume of the ring counter is less than that of the central counter. 

RECOIL-SPECTRUM CONSIDERATIONS 

It should be much easier to draw some definite conclusions about 
the suitability of the tubers multiplication factor for obtaining usable 
pulse heights if a monoergic neutron source were used, preferably one 
giving a recoil alpha in the energy range of interest, 770 kev to 1.80 
Mev. The counting results from a monoergic source would be much 
easier to analyze. 

This is illustrated as follows. For a monoergic source the 

neutron spectrum, N(E^) versus E^, is ideally a spike and N(E^) is 

the total number of neutrons; the differential pulse-height spectrum 

due to the recoils, n(W) versus W, where W is pulse height, is a 

rectangle; and the integral bias recoil spectrum, N(W) versus W, is 

ideally a constant- slope line. These are sketched in Fig. 7. These 

shapes of the differential and integral bias curves are justified as 

follows. The pulse height depends on the energy of the recoil nucleus; 

the counter and amplifier convert that energy into pulse height. The 

recoil nucleus, as shown in Chapter VIII, can have a maximum energy 

2 

of K E , where K =4 A/(1 + A) , if we use> mass numbers for the 
a n^ a ' 

masses. Let W be the energy of recoil nucleus (converted into volts 
of pulse height by the counter system). Then W can vary from 0 to 
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Fig. 6o Integral bias curves for ring counter at various 
counter voltages, 

A. 10Z5 volts, 

B. 975 volts, 

C. 9Z5 volts, 

and using neutron source in arrangement of Fig. 4. 
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7. Ideal curves for monoenergetic neutron spectrum. 

A. Neutron spectrum. 

B. Differential bias recoil spectrum. 

C. Integral bias recoil spectrum. 
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and the number of pulse heights from recoils from the monoergic 
source are spread out over this range of W values. The question is 
whether or not the distribution is uniform over the range, giving a 
rectangle. At the values- of of interest here, and for light nuclei, 
it is accurate to assume that neutron scattering is mainly S-wave and 
spherically symmetric in the center- of-mass system of coordinates. 

We wish then to express the recoil energy in terms of E^ and the 
angle of scattering in center-of-mass coordinates (both energies in 
laboratory- system coordinates). Using conservation of energy and 
an equation developed** ^ for relating the velocities of the neutrons 
in lab coordinates before and after scattering with 0, the angle of 
scattering in c.m. coordinates, we can show W = ^‘K^E^ (1-cos 0); 
then dW = K E sin 0 d0. But in spherical coordinates, we have 
d0 = df^/ZTT sin 0 , where dQ = element of solid angle, therefore 
dW = (dn/4TT). We can say Prob(W) dW = {<j(6)/q dfi , 

where 

a(0) = differential scattering cross section (c.m. ), 

a ^ = integral scattering cross section, 

Prob('W) dW = probability for recoil nucleus to acquire an energy 
between W and W + dW , 

(cr ( 0 )/c 7 g)df2 = probability that neutron will be scattered through 
angle 0 into dQ . 

It follows: Prob (W) = cr(G)/(J ^ * 4tt/K^ E^ . F or an isotropic 

scattering distribution (a ( 9)/ ^ ^ ^ 1/4tt) and monoergic neutrons, 

Prob(W) is thus a constant, 1/K E , and n(W) is a rectangle of 

a n 

constant value proportional to N(E^), The integral bias distribution 
is an integral of the rectangle, integrating from maximum W to zero, 
giving the triangular shape. 

For a non-mcnoergic neutron spectrum, N(E^)- -instead of being 
the total number of neutrons -- represents the number of neutrons per 
energy interval, and we must integrate over the spectrum. We can 
say 
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n(W) K, r Prob (W^ N (E » ct E | dE K, f L/K E N'E ' criR )dE„ 
1 J n s n 1 j ' a n n 



where K. is an arbitrary constant and a (E ) is ^he elastic scattering: 

cross section as a function of neutron energv„ Thus for a complex 

neutron spectrum the shapes of the differential and integral bias 

curves are very difficult to predict. For perform'mg the integration 

in analytical form an expression for N E ! and cr E Imust be used 
^ ^ n s n*^ 

Generally, for most sources N'^E ]) goes from zero at E 0 to a 
^ n ^ 



n 



maximum and back to zero at E - E (maxK The intervening 
structure may be quite complicated. Let us assume J E ) varies as 

I i 8 ] ^ ^ 

This IS approximately true'- ■' for scattering on protons 
an approximate relation being (.or ) ^ 



1/\JE 



n 



for E between 0.2 and 5 Mev 
n 



4 . 5 / n/ E for E in Mev and a in barns„ It is less true for heUum, 



n 



which has a peak of about 6.b barns at i.iS Mev and then tails off 
Then if, for illustrative purposes, we assume 

3/2 



N(E ]) - E 
n^ n 



sin E 

n 



where E goes from zero to E (max) 



n ^ n 

a spectrum shape that can be readily integrated 



n arbitrary units we ha e 



E 

n 

n(W) = K, f l/(K E ) E sih E H ~ ^ dE 

'' ''an'n nn n 

E =W/K 
n ' a 

K, 

= -j— (1 + cos (w/k^):;, 

a 



where W goes from 0 to K it. To save writing W/K as the variable 

a a 

on the next integration, let us change W from energy to arbitrary 

voltage units where W goes from 0 to it. Since we have divided the 

W scale by K we must multiply the density d^strrbutibn by K . The 
a a 

differential distribution becomes 

n(Wi = Kj r cos 
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The integral bias distribution is obtained by integrating this from 
the right: ^ 

N(W) - - W) dW = (tt - W - sin W). 

The value of this distribution at W = 0 is the total number of scattering 
events, This should equal the number of neutrons times the 

probability of scattering or 

'TT TT 

K, ( N(E )a ‘ (E )dE =K, f E sin E dE = ttK, . 
lj'n*s'n'n Ijn nn 1 

0 0 

The three distributions are sketched in Fig. 8. The significant 

trend is the way the integral bias distribution rises rapidly near zero 

pulse height, even with this assumed neutron spectrum weighted 

heavily toward the high-energy end of the spectrum. With a complex 

neutron- energy spectrum--as obtained, for example, from (a, n) 
f 37l 

sources --the shape of the integral bias curve and the manner in 
which it tends to the zero pulse-height counting rate becomes very 
difficult to predict. 
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Fig. 8. Assumed neutron spectrum and recoil spectrum 
for assumed conditions. 

A. Neutron spectrum. 

B. Differential bias recoil spectrum. 

C. Integral bias recoil spectrum. 
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CHAPTER VI 



PARTICLE RANGE AND WALL EFFECT 

Kaowing the range of an alpha particle or proton in air, we can 

compute the range in helium if we know the relative atomic stopping 

power of helium compared with air. The range-energy curves to be 

used for alpha particles and protons in air at 15^ C and 760 mm 

f 3 8 1 

pressure are those of Bethe, ^ ^ which, on the basis of experiment, 

corrected for low energies the original curves by Livingston and 

[ 39 ] 

Bethe , The original work was a theoretical treatment based on the 
Born approximation but also modified by experimental data then 
available. 

The basic equation giving energy loss per cm of path defines 
stopping power: 

- dE/dR = (4tt e^z^mv^) NB, 

where 

E = energy of incident particle, 

R = path length, 
e = electron charge, 

z = atomic number of incident particle, 
m - electron mass, 

V = velocity of incident particle, 

N = number of atoms per cc of stopping material, 

B = stopping number = Z In (2m v^/l ), where 

a V 

Z = atomic number of stopping material, 

= K Z = average excitation potential of atom of stopping 
material, 

K = Bloch constant (which is not actually a constant but varies 

somewhat for different atoms and is usually determined 
experimentally because it is difficult to calculate 
theoretically). 

The relative stopping power, compared with air, is defined: 
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s B/Bq - Z/Zq[(1!i 2mv^ - In I^^)/(ln 2mv^ - In 1^)] ^ 



where the zero subscripts refer to air. Here s is seen to be a function 

of incident-particle velocity, and for different atoms increases more 

slowly than Z alone. We find it accurate enough for the purpose at 

hand to use this equation, although many modifications to the theory 

can be made in computing B, as done by Hirschfelder and Magee^^^^ 

for several different substances for protons. 

We take experimentally determined values of air from 

Reference 39- ~ 80.5 ev and Z^ 7.Z2. For helium I = 44 ev^"^^^. 

0 0 av 

The calculation shows, for 1-Mev protons, ^ 0.328. To convert 

the range in cm in air at NTP (15° C, 760 mm) to that in helium under 
the same conditions we note that stopping power is proportional to 
N and B, therefore the range in cm must be inversely proportional to 
these quantities. The range in helium is 

R = Rq • Nq/N • Sq/s^^ , 

where the zero subscripts again refer to air and Sq - 1. Now, Nq/N 

is approximately 2, since air is mostly diatomic. The calculation 

yields R = 6.13 Rq. This should be fairly accurate for helium-3 or 

helium-4. For a 1-Mev proton, Rq is 2.28 cm, therefore R is 14 0 cm. 

The range is inversely proportional to gas pressure, since density is 

proportional to pressure. Thus, at Z.5 atmos this proton range in 

helium becomes 5.60 cm, and at 5 atmos, 2,80 cm. 

To find the range of an ion of other isotopes of hydrogen (given 

the relation for protons) we note that of all terms in the stopping 

Z 

power formula only v is different for a given energy. Thus, a 

triton of energy E experiences the same stopping power as a proton 

of energy l/3 E; but it goes three times as far since it has three times 

the energy. The maximum triton energy with which we will be con- 
3 

cerned in the He disintegration is about 1 Mev, and we figure its 

range to be 8.54 cm in 1 atmos of helium. Since this is less than a 

3 

representative proton range from a He disintegration and represents 
an upper bound on the triton range, the order of magnitude of the 
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wall-effect difficulties can best be estimated by considering the 

[ 7 ] 

proton range. As ointed out by Batchelor in explaining his computer 

program for wall-effect computation, the correlation between proton 

and triton tracks plus a cylindrical geometry makes exact calculations 

3 

impracticable; also the angular distribution of the He (n, p)T reaction 
IS not very well known. 

In integral bias counting involving the (n, p) reaction, detection 
efficiency is defined as the counting rate divided by rate of emission 
of protons, When the total range of particles is small compared with 
linear dimensions of a counter and distribution is isotropic, a value 
for detection efficiency often given is 

F(E) = 1 - R^/2 b - L, 

where we have an active cylindrical volume of radius b and length L, 
and is that portion of total range of protons necessary to produce 
a pulse above the bias energy E of the discriminator ^ . The terms 

in the expression for F(E) reducing it from unity represent the 
magnitude of the wall effect and indicate the loss at the lateral bound- 
ary and end face, respectively. 

The expression for F(E) can be shown to be plausible as follows: 
Consider the fraction of the volume of a cylinder of radius b within a 
distance R of the side walls. This fraction is (2 bR - R )/b . 

Consider charged particles originating with equal probability along a 
radial line with ends defined by the outer volume. If any direction of 
emission is equally probable, one-half of the particles will be emitted 
with a radial component toward the center and not intercept the wall. 
For that half of the particles having a positive radial component, 
those originating on the inner edge of the fractional volume have 
practically zero probability of intercepting the walls (since only one 
direction of emission allows the particles to reach the wall), while 
those originating at the wall have 100% probability of reaching the wall. 
Assuming that this increase in probability is linear as we move out- 
ward toward the wall, we estimate that one -fourth of all the particles 
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originating in the fractional volume reach the wall. Then the wall 
effect is l/4 time the fraction of the volume, or 

R/2b - (R/2b)^ . 

For R 2b, this is approximately R/Zb. Similarly, the fractional 

volume within a distance R of the ends of a cylinder is 2R/L, and-- 
again taking l/4 of this value-- we get R/ZL as the wall effect due to 
the ends. 

In our case the condition of large dimensions is not well ful- 
filled; also we assumed that the ionization from the particle other 
than a proton was confined to a very small region, which is not so 

3 

true for a He reaction as for one involving heavier nuclei. In 
differential detection F(E) is ideally the difference between two step 
functions with different energy thresholds. Therefore, the deviation 
from unity is a rough indication of the magnitude of wall effect on 
energy determination for a proton of given energy. For a 1-Mev 
proton in 2.5 atmos of helium, the geometric factor is 

G(E = 1 Mev) = I - R/Zb - R/ZL, 
and with Zb “ 7.62 cm, ZL = 38.9 cm, R = 5.60 cm, we have 
G(E)“0.1Z or the wall effect is 88%. Doubling the pressure would 
reduce the wall effect to 44%, and at 10 atmospheres to 22%. 

For counting neutrons in ordinary helium we are concerned 
with the range of the alpha recoil. The maximum energy of the alpha 
particle is 16/Z5 E^, as shown in Chapter VIII. If we consider; the 
average energy of the Pu-Be neutron spectrum to be 4 Mev, this 
corresponds to a maximum alpha energy of 2.56 Mev, or an average 
E^ of about half this amount. Figuring the range of this alpha 
particle in 2.5 atmospheres of helium by our previous method, we 
get 1.38 cm. The wall effect is 22% for this range. Experimentally, 
a coincidence counting rate of about 16% of the total counting rate in 
the main counter was measured. This is a very rough check, but it 
at least appears that the ring counter is capable of reducing the wall 
effect. 
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ELECTRON RANGE 



To obtain the range and energy loss of electrons in estimating the 

effect of tritium decays or background gamma counts, one can take the 

[41 1 

nearest tabulated value*- of stopping power for helium as a function 
of electron energy, given in units of Mev per g/cm . For example, at 
10 kev, we have 

- dE/dx = 22.5 Mev/g/cm^ c 

If we take A x as the average path in the counter, 1.59 cm, times the 

-4 / 2 

density of helium at 2.5 atmospheres, Ax = 7.53 • 10 g/cm . Then 
AE IS greater than 16 9 kev (the product of - dE/dx ' Ax), because 
the stopping power increases with decreasing energy. For higher- 
energy electrons it is fairly accurate to use Feather’s Rule or 
Flammer sfeld’s formula to find x and then multiply Ax/x times 

the energy of the electron. 
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CHAPTER VII 



NEUTRON GENERATOR 

•T"or testing purposes one desires a monoenergetic neutron 

source. The complex spectrum of (a, n) sources has been mentioned; 

the other class of neutron sources using radioactive nuclei are the 

photoneutron or {y, n) sources. These use deuterium or beryllium 

because other nuclei have thresholds above 6 Mev. To get high enough 

intensity the gamma source is surrounded so that there is an intrinsic 

energy spread due to different angles between neutron emission and 

incident gamma ray. This is small; for example, it is 3.37.% for a 
24 9 fzOl 

Na + Be source*- . However, this is not the main cause of 
energy spread. Owing to the large quantity of beryllium there is 
considerable neutron scattering in the source, and gamma-ray loss 
of energy by Compton scattering. 

It has been found possible to construct neutron generators 

3 

utilizing the deuter on- deuteron reaction, D(d, n)He , with useful 
neutron intensities without elaborate high-voltage equipment; for 

[421 

example, the generator of Zinn and Seely was built in 1937. 
Voltages on the order of 100 kv can give useful yields even though the 
reaction is not very efficient at this energy. With this low voltage we 
also have the advantage of almost monoenergetic and gamma-free 
neutrons. For a thick target the deuter on- energy at the time of the 
reaction can vary from a maximum to zero. In the generator used 
in this experiment the maximum deuteron energy is 120 kev. Thus, 
using equations developed in Chapter VIII and Q “ 3.28 Mev, we find 
that in the forward direction the neutron energy can vary from 2.46 
to 2.91 Mev. At a laboratory- sytetem angle of 90 degrees the energy 
spread is from 2.46 to 2.49 Mev, or little more than 1%. The energies 
are weighted toward higher values because yield increase s rapidly 
with deuteron energy. This neutron- energy range is convenient for 
checking the operation of the counter because the maximum recoil- 
alpha energy is about the same as the maximum (E^ + Q) energy we 
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wish to measure in helium-3, and about half the recoils are above the 
lowest sensitivity needed. 

PULSED- NEUTRON GENERATOR 

The neutron generator that is being put into operation provides a 

pulsed source of neutrons and is similar to the one described by 

f43l 

Ruby except for the method of mounting the accelerater system and 

the vacuum equipment. The accelerator system is described in 

Reference 43; briefly, the deuterons are generated in a pulsed Phillips 

Ion Gauge type of ion source supplied with deuterium gas through a 

palladium leak and accelerated to impinge upon a circular target of 

l-l/8-inch diameter. The targets are prepared by occluding deuterium 

in titanium. This type of target allows simpler equipment than one 

with heavy ice, by requiring no special cooling. It is planned to try a 

deuterated paraffin target to see if the yield is increased. The 

12 13 

objection of non-monoenergetic neutrons due to the C (d, n)N re- 
action is not pertinent here because we are below the reaction 
threshold. The generator was designed so that adjustment of param- 
eters would not be critical; the yield is largely determined by the 
accelerating voltage and the condition of the target. (The reaction 
cross section decreases with deuteron energy, and the targets have a 
tendency to form a carbonous surface after a few hours* use. ) 

The pulsed operation is particularly convenient for experiments 
that require a reference time, such as in neutron-diffusion experiments; 
single pulsing has been used for cloud chamber experiments. In normal 
operation the ion source is pulsed at 60 cycles per second with 200- 

psec pulses. In this manner we can get an integrated yield, Q , of 

6 ^ 

more than 10 neutrons per sec. This is a duty cycle of 0.012, or 

during a pulse the instantaneous neutron yield is Q^/O.012 neutrons 

per second. For a cloud chamber experiment we used an external 

oscillator option of 400 cycles per second and increased the width of 

the pulse by an adjustment in the arc pulser to 500 psec. This is a 

duty cycle of 0.2. This mode of operation was gated on, when it was 
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desired to take a picture, for 100 millisecoads, giving a total pulse-on 
time of ZO msec per gate. Then the neutron yield per gate should be 
Qg/O.OlZ times ZO msec, or more than 10^ neutrons in each burst, 
assuming the same effectiveness of the ion source when pulsed at 400 
cycles per sec. as at 60 cycles. It would be more realistic to calibrate 
at 400-cycles steady operation, but this increases the risk of damage 
to the ion source. 



NEUTRON YIELD 

In order to measure the integrated yield one should determine 
if the anisotropy of the yield is significant. It has been found that in 
the energy region below 0.5 Mev, the neutron yield as a function of 



angle in center-of-mass coordinates can be represented 



[44] 



by 



N( e) = A (1 + B cos e), 

where A is a constant and B is a number that increases with deuteron 
energy. Thus, the yield is smallest at 6 “ 90 degrees. The quantity 

[451 

B has been determined experimentally ^ down to 500 kev, and the 
curve extrapolates fairly well into a theoretical determination^^^ ^ at 

3 

low energies based on the D(d, p) H reaction. It is believed the 
angular asymmetry is the same for the companion reactions at 
corresponding voltages. We estimate B as 0.6 at IZO kev. This 
means the neutron flux at 0 = 0 should be (1 + B) times the neutron 
flux at 6 - 90 degrees. Or, performing an integration, we find that 
if is determined by measurements at 6 = 0, we should divide this 
by (1 -KB/Z+B)) for the true yield or 1.Z3 for B “ 0.6. This is a high 
estimate of the anisotropy, because not all the deuterons have the 
maximum energy. With these small deuteron energies, we expect 
that the labor atory- system angles will not be much different from the 
center-of-mass angles. We need not outline the calculations here 
(a description of the two coordinate systems and some convenient 
equations may be found in Reference 47), but for 0-90 degrees 
we find the lab angle ® “ 85.6 degrees for IZO kev, and for 

degrees, 0 = 94.4 degrees. 
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The generator yield has been monitored with a proportional 

counter surrounded by 2.5 inches of lucite moderator and a thin sheet 

of cadmium. Care must be taken in its calibration since this 

arrangement does not give equal sensitivity over a wide range of 

neutron energies. For a particular counter we found a sensitivity of 

2 

0.29 counts per neutron per cm when a Pu-Be source was used, and 

2 

0.37 counts per neutron per cm when a mock fission source with a 
smaller energy range of neutrons was used. However, it can be 
calibrated on the D-d neutrons by comparison with a calibrated 
Hanson and McKibben type^^^^ long counter. The lucite embedded 
counter is easier to place and adds less scattering objects to the 
experimental arrangement. 



- 49 - 



ft lUi** X — ■. » trni^Jitr-mt^tmmrnm »m0f 

*»*m 

• ftni MIM ttt t» iwfcit <>y ft— 1» ■■ w^f npi’wwM I* 

ftrtft' ftftfti « I r «0 fnVlrtftW^* » 1 #»-- " 11*11 1 »ft»#M»»ftft 

•4 ^ I Uaivi I §g|' * 90m ftHMHMft 

IftHU jm*my»^P -• • ^ ^ 

• i1l» nr«anii» ■»•••• j«m| # ■— ••■MIC ifvO 

• ItlMfl' mty laftt iiXiAlMlI.Yi'l fW*| ‘'^ftftM 

dHMvnrrN^ • .ifrtt ^ftvftmrii . ^ -9 Ifttift *11)1111^ 

iVjtfirM«n<«l|>ll«uPl «fi» 111®' 4 «W0I If' •- 

ig^ Mg|)liMl Hi^ni • I •• # ««MHi4 

ivrftvwi iw fft I 



CHAPTER VIII 



NUCLEAR REACTION DYNAMICS 
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Consider the general case Sketched above of a particle of mass 
and velocity colliding with a target nucleus of mass assumed 
at rest. After the reaction the resulting particle of mass flies off 
with velocity v^ at angle The nucleus of mass has velocity 

v^ at angle We will use nonrelativistic dynamics, mass numbers 

for masses, and laboratory coordinates. 

First, write two equations for conservation of momentum for 
components parallel and perpendicular to Vj^: 



^ 1^1 ” ^ 4^4 ^ ^ 3^3 



(VIII- 1) 



M 3 V 3 sin 0 j^ ^ ^ 4^4 



(VIII- 2) 



By transposing the last term of Eq. ( 1 ) to the other side, 
squaring both equations and adding, we can eliminate Then 

replacing velocities by use of E == M v for energy and rearranging, 
we obtain 



M 



1 




E 



4 M 



4 
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We can write another equation for 
Q of the reaction. Having assumed E^ = 0, 

E^ = Q + Ej-E3. 



from the definition of the 
we write 



For the He (n, p)T reaction, using Q = 0.77 Mev and similar 
symbolism, we write two equations for the triton energy: 



Et = ^ En + ^ Ep - I (E^ Ep)2 cos 



e 



L ’ 



E^ = 0.77 + E - E 
1 n p 



With these two equations we can solve for the proton and 

triton energies for given values of E and 0^ . 

n X-/ 

F or = 0, the proton receives 75% of the Q energy and the 
triton receives the remaining 25%. That is, 

E = 0.577 Mev, 

P 

E^ = 0. 193 Mev . 

For E^ = 1 Mev, a tabulation is given below. 



(degrees ) 


Prot 


on 


Triton 


E 

P 

(Mev) 


% of 
(E^ 4 Q) 


(Mev) 


% of 

(E^ + Q) 


0 


1.74 


98.3 


0.03 


1.7 


45 


1. 52 


85.8 


0.25 


14.2 


90 


1.08 


61.0 


0.69 


39.0 


135 


0.77 


43.5 


1.00 


56.5 


180 


0.67 


37.8 


1.10 


62.2 



It is often desirable to have a convenient expression for the 
energy of the ejected particle at 0 = 0 and at 90 degrees as a function 

of the bombarding particle energy. For the forward direction the two 
equations become 
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E, - (E, E, M, MJ2 , 



1 



1 



E 4 = Q + E^- E 3 . 



A solution for E yields 



E3 = A3/2A3 + i [(A2/A3)^/ - ( 4 A^/A 3 )] 2 . 



where 



A^ = M^Q + E^ (M^ - M^) , 

A, = 2A^ (M, + M^) + 4 E, Ml M, , 



1 ^”1 



A 3 = (M 3 + M^)^ . 



For the He (n, p) T reaction this gives 



n 



For the D(d, n) He reaction, 



^n = |0" i ^d<‘ * ^13 + 



]M. 



For 0j^ = 90 degrees, the two equations are 



- M M. 

^4 = M7 + MT ^3 • 

4 4 



£4 = Q + E, - E 3 



Solving, we obtain 



^3 = 



M4/(M3 + M^) 



Q + [(M^-Mj)/(M^ + M 3 )] E^ 
For the two reactions above, at 90 degrees. 
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3 1 

E = 4 Q + 1 E , . 

n 4 4 d 



For the case of elastic scattering we set Q = 0 and = M^, 
therefore the two equations become 



M 

^4 = m' 
4 


Ml 

M 4 


E4 = Ei-E3. 





2 

M, 



(E^ E^ M^)2 cos , 



We are interested in the energy given to the recoil nucleus, 
therefore we eliminate from the first equation by using the second. 
Then we set 0^^ = 180 degrees for maximxim E^ corresponding to a 
backward recoil. Solving for E^(max), we have 

4M M 

E (max) = ^ E, ; 

For neutron scattering in helium- 3, 

EHe3 (max) - | ; 

and for neutron scattering in helium-4, 

E (max) = ^ 

a ' 25 n 
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